X(3915) and X(4350) as new members in P-wave charmonium family 
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The analysis of the mass spectrum and the calculation of the strong decay of P-wave charmonium states 
strongly support to explain the newly observed X(3915) and X(4350) as new members in P-wave charmonium 
family, i.e., x' c q f° r 3TC3915) and^-" 2 for X(4350). Under the P-wave charmonium assignment to X(3915) and 
A"(4350), the J pc quantum numbers of X(3915) and A"(4350) must be ++ and 2 ++ respectively, which provide the 
important criterion to test P-wave charmonium explanation for X(39 1 5) and X(4350) proposed by this letter. The 
decay behavior of the remaining two P-wave charmonium states with the second radial excitation is predicted, 
and experimental search for them is suggested. 

PACS numbers: 14.40.Pq, 13.25.Gv, 12.38.Lg 



Recently two new charmonium-like states X(3915) and 
X(4350) were released by the Belle Collaboration in the yy 
fusion process Q H. X(3915) is observed by the invari- 
ant mass spectrum of J/iJ/u) in yy — > J/if/ai channel. The 
mass and width of X(3915) are M = 3915 + 3(stat) + 2(sys) 
MeV and r = 17 + 10(stat) ± 3(sys) MeV (J. X(4350), a 
new charmonium-like state found in the invariant mass spec- 
trum of J/iff(p, is of m = 4350+^(stat) ± 0.7(sys) MeV and 
F = 13.3^ 7 j 9 (stat) + 4.1(sys) MeV J!. 

Until now, Belle experiment has reported three 
charmonium-like states via the yy fusion. Besides X(3915) 
and X(4150), Z(3930) is a charmonium-like state observed 
in yy — > DD, which is of mass m = 3929 + (stat)5 ± 2(syst) 
MeV and r = 29 + 10(stat) ± 2(syst) MeV Jl. The angular 
distribution in the yy center of mass frame shows J PC = 2 ++ , 
which indicates that Z(3930) is a good candidate of x' c2 , i.e., 
a charmonium with « 2s+1 7l = 2 3 f2 Jl- 

The observations of X{39\5) and X(4350) ont only make 
the spectroscopy of charmonium-like state observed by the 
yy fusion process become abundant, but also help us further 
reveal the underlying structure of charmonium-like states ob- 
served by the yy fusion. Before illustrating the underlying 
structure of X{39\5) and X(4350), one first gives a brief re- 
view of the established P-wave charmonium states or their 
possible candidates by Fig. Q] where /i e (3525) (a P-wave state 
with spin 0) is not listed. Three P-wave states without the ra- 
diative excitation are^ c o(3415),^ e i(3510) and^ C 2(3556) 0]. 
For the first radial excitation of P-wave charmonium, the can- 
didate for ++ state x' c0 is still absent while X(3872) JlH3l and 
Z(3930) Jl 

can be recommended as 1 ++ state x' c \ aR d 2 ++ 
state x' c2 i respectively. 

Usually the yy fusion process provides a good environment 
to create charmonium by yy fusion into a pair of cc. Although 
Z(3915) seems to be explained as an exotic state indicated in 
Ref. JH], in this letter we adopt another point of view to ex- 
plore whether the newly observed X(3915) can fill in the blank 



of the remaining P-wave charmonium with the first radial ex- 
citation shown in Fig. Q] The mass of X(3915) is consistent 
with the result of the potential model, which once predicted 
the mass of the first radial excitation x' M is around 3916 MeV 
according to Godfrey-Isgur relativized potential model 
Since there exists the vanishing coupling of X(3915) - DD* 
and the week interaction of Z(3930)-DD* while X(3872) with 
J PC = 1 ++ interacting with DD* via S-wave is very strong, the 
coupled channel effect on bare;^ and^-, is weaker than that 
on bare x' c \ JH 0]> which explains why the mass difference 
between X(3915) and Z(3930) is smaller than that between 
X(3915)and X(3872). 
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FIG. 1: (Color online.) The established P-wave state charmonium 
states without the radial excitation 01 and the candidate for the first 
radial excitation of P-wave state. A comparison between newly ob- 
served X(4350) and the candidate of P-wave charmonium states is 
given. 



By the comparison of X(4350) with the existed P-wave 
states, one notices that the mass difference between X(4350) 
and Z(3930) is about 420 MeV, which is similar to that be- 
tween Z(3930) and x C 2(lP). The regularity of the mass gaps 
existing xa(lP), Z(3930) and X(4350) is consistent with the 
estimate from Resonance Spectrum Expansion (RSE) model, 
which indicates the mass gap between the states with the ra- 
dial quantum numbers « and « + 1 is 380 MeV JIolfTlll . Thus, 
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one further proposes that X(4350) is as the second radial ex- 
citation of P-wave charmomium state. We also notice the pre- 
diction of the mass of x" 2 in Ref - H b y Godfrey-Is gur rela- 
tivized potential model, which is about 4337 MeV. This value 
is consistent with the mass of X(4350). 
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TABLE I: The allowed open-charm strong decays of x' c0 an <ix"j (J - 
0, 1,2). Here, we take 4350 MeV as the upper limit of the mass 
of x"j- £>i(2420) is the 1 + state in the T = (l + ,2 + ) doublet while 
D i (2430) is the 1 + state in the 5 = (0 + , 1 + ) doublet since as indicated 
in Ref. 03]. 



The decay modes of the created charmonium state from 
the yy fusion include open-charm and hidden-charm de- 
cays, which are the observed decay channel of Z(3930) and 
X(3915)/X(4350), respectively. For testing the proposal for 
the structure of X(3915) and X(4350), in the following, we 
further study open-charm decay of the radial excited P-wave 
charmonium x' c q and y 1 ', ( J - 0, 1 , 2) by the Quark Pair 
Creation (QPC) model lll2]4T4il . which is a successful phe- 
nomenological model to calculate Okubo-Zweig-Iizuka (OZI) 
allowed strong decays of hadron. 

The allowed decay modes of y' M and x", are presented in 
Table Q] The QPC model provide us an effective approach to 
study the two-body strong decays of the radial excited P-wave 
charmonium x' c q an< ^Xcj- 



The transition matrix element of the process A(c(l)c(2)) — » 
B(c( l)q(3)) + C(c(2)q(4)) in the center of mass frame 
of charmonium A is written as {BC\T\A) = <5 3 (K# + 
K c )M Mj a Mj b Mj c(K), where the transition operator T in the 
QPC model reads as 

T = -3y^(lm;l - m\Q 0> f dk 3 dk 4 5 3 (k 3 + k 4 ) 

xJ/ lB! (^i) X \U # < 4*3) b\/M , 
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TABLE II: The partial wave amplitude of the open-charm decay of 
x' c0 andx' c 'j- Here, two 1 + charmed mesons in S and T doublets are 
the mixture of two basis 1 1 Pi and l 3 Pi, which indicates |l + (5)> = 
coseil'P,) + sineil^,) and |l + (r)> = - sin6(|l 1 J Pi) +cos6>]l 3 P,) 
03l . One takes f = 1/ V3 obtained from the calculation of the 
flavor matrix element. In heavy quark limit, one usually takes the 
mixing angle 8 = -54.7°. According to the partial wave amplitude, 
the partial decay width is expressed T = n 2 \M jl \ 03-03]- 



where i and j denote the SU(3) color indices of the created 
quark and anti-quark from the vacuum. ^ = (uu + dd + 

55)/ V3 and w 34 = d a3a Jy*3(a = 1,2,3) means flavor and 
color singlets, respectively. X^\- m i s a triplet state of spin. 
J/f,„(k) = |k| f YfmiOk, 4>k) is the (th solid harmonic polynomial. 
Dimensionless constant y is the strength of the quark pair cre- 
ation from the vacuum and can be extracted by fitting the data. 
In this letter, y = 6.3 lfl5ll . The strength of ss creation satisfies 
y s = 7/ V3 111. The helicity amplitude M Mj * Mj s m 'c (K) is 
extracted by the transition matrix element, which is related to 
the partial wave amplitude by iflrjl 



M JL (A^BC) = ^ L+1 V (L0JMj a \J a M Ja ) 

x{J b M Jb J c M Jc I JM Ja )M M] « Mj e M 'c (K), 
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where J = J B + Jc and J A + Jp = Jb + Jc + L- A detailed 
review of the QPC model was illustrated in Ref. Il7tl . The par- 
tial wave amplitude corresponding to the open-charm decays 
shown in Table U is presented in Table [TT] The concrete ex- 
pressions of Ojj, Qij and Vij,k are extracted from the spatial 
integral I^'m (K), which describes the overlap of the ini- 
tial meson (A) and the created pair with the two final mesons 
(B and C). Here, the harmonic oscillator (HO) wave function 
fnrfmCk) = fl nr ((R, k)J/„.f„,(k) is involved in the calculation of 
the spatial integral. Parameter R in the HO wave function is 
obtained by reproducing the realistic root mean square (RMS) 
radius by solving the schrodinger equation with the linear po- 
tential 111511 . The relevant mass values are taken in PDG 14[]. 
We take R = 1.52, 1.41, 1.85, 1.69, 1.85, 1.75, 2.00, 2.00 
GeV 1 corresponding to D, D s , D\ D*, D*(2400), D* (2317), 
Di(2430) and Dj(2420), respectively. Besides, other parame- 
ter inputs include m c = 1.6 GeV, m u — mj = 0.22 GeV and 
m s = 0.419 GeV. 

In the right diagram of Figs. [2] one presents the total width 
of/f'n w ^ m tne variation of R. The node effect from the wave 
function of higher radial excited states results in the decay 
width calculated by the QPC model being dependent on the 
R value. When taking R = 1.80 ~ 1.99 GeV 1 , the obtained 
total open-charm decay width of x'm falls in the range of total 
width of X(3915) released by Belle JT|] (the calculation result 
with the typical value R = 1.92 GeV 1 corresponds to the 
central value of the width of X(3915)). As &x' C 2 charmonium 
state JH, Z(3930) can be as a realistic test of the reasonability 
of the range of R for;^ . Our theoretical result of the open- 
charm decay of x^ dependent on R and the comparison of the 
calculation result with Belle data of Z(3930) yj], which are 
shown in the left diagram of Figs. [2] indicate that the upper 
limit of R for;^ 2 is very close to the lower limit of R for^ as 
marked by the red arrows in Fig. |2j which further shows the 
reliability of investigating;^ open-charm decay in the range 
of R = 1.80 ~ 1.99 GeV- 1 . Thus, explaining X(39 15) as ax' c0 
charmonium is tested through the open-charm decay of x' c(y 

Considering an exotic D*D* molecule explanation to 
X(3915) suggested in Ref. |8y, one proposes the experimental 
study of the open-charm decay DD to be as a best way to dis- 
tinguish between the exotic and the conventional states for the 
controversial X(3915) since the DD decay of X(3915) under 
the D*D* molecule assignment occurs via hadronic loop effect 
as indicated in Ref. 1181 . which results in the decay width of 
X(3915) — > DD under the D*D* molecule assignment being 
far smaller than that under x' c0 explanation to X(3915). 



The results of the open-charm decays of x"j (J = 0, 1,2) 
are presented in Fig. [3] via scanning the parameter space 
R = 1.8 ~ 3.0 GeV 1 , which is due to the radius R of x' c 'j 
being fatter than that ofx' cJ - We need to emphasize that the yy 
fusion process determines the most possible quantum number 
of X(4350) to be ++ or 2 ++ , which makes us choose X(4350) 
as the candidate of x"q or x"i an d fully exclude x'U assign- 
ment to X(4350) 10]. The open-charm decays of x'Jn and;^', 2 , 
two candidates of 7(4350), display different behaviors as il- 
lustrated in the left and right diagrams of Fig. [3] The total 
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FIG. 2: (Color online.) The dependence of the decay width of 
Z(3930) and X(3915) on R under X ' cl and X ' c0 assignment for Z(3930) 
and X(3915), respectively. Here, red dash line with grey band denote 
the central value for the error of total width of X(3915) and Z(3930) 
measured by Belle yl§]. The green band denotes the region of R re- 
sulting in the theoretical values consistent with Belle data. The solid 
lines with blue error bands are our calculation result. 



open-charm decay of^" 2 with R = 1.9 ~ 2.3 GeV is well 
consistent with Belle data J2l as shown in the right diagram 
of Fig. [3] while the total open-charm decay of x" Q is far away 
from the Belle data H, which shows that we can fully exclude 
x' c ' explanation for X(4350) and finally establish X(4350) as 
a good candidate of x' c ' 2 - Our numerical result demonstrates 
that rr/04" (DD* + h.c.,D s D* + h.c, D*D*, D*D*) are 
the dominant decay channels for X(4350). 
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FIG. 3: (Color online.) The variation of the decay width of x"j (J - 
0, 1, 2) with R value. Here, we set the upper limit of the masses of P- 
wave states with the second radial excitation as 4.35 GeV. The yellow 
dash line and shaded grey band shown in diagram of %" 2 denote the 
central value for the error of total width of X(3915) measured by 
Belle d. 
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Meanwhile, the predicted properties of the remaining two 
P-wave charmonium states with the second radial excitation 
can be as the guidance of experimental search for^',' () and^" r 
The predicted total decay width of %" is around 82 ~ 110 
MeV corresponding to R = 2.0 ~ 3.0 GeV which is not 
strongly dependent on the R values. As the dominant de- 
cay mode of 0~0~ channel including DD and D+D~ is 
a golden channel to find^" . For^" p its total decay width 
is of large span from 47 MeV to 140 MeV corresponding to 
R = 2.0 ~ 3.0 GeV 1 . Among its partial decay channels, 
0~1~ channel is always the main decay channel of%" 1 under 
taking different values of R, which indicates that DD* + h.c. 
and D + S D*~ + h.c. can be as the suggested decay channel of 
searching for^" r The detail of the open-charm decay behav- 
iors ofx' c 'j are listed in Fig. [3] 

In summary, the newly observed X(3915) and X(4350) are 
firstly explained as an d x" c i extremely well respectively 
by analyzing the mass spectrum of P-wave charmonium fam- 
ily as well as by calculating the open-charm strong decay of 
X(3915) and X(4350), which are consistent with the existed 
experimental findings. Just because of our explanations to 
X(3915) and X(4350), the spectroscopy of P-wave charmo- 
nium becomes abundant. Under the assignment of P-wave 
charmonium to X(3915) and X(4350), the J PC quantum num- 
bers of X(3915) and X(4350) will be definite, i.e. J PC = ++ 
for X(3915) and J PC = 2 ++ for X(4350), which provides the 



powerful criterion to test the P-wave charmonium assignment 
for X(3915) and X(4350) since the experimental analysis of 
the angular distribution of X(3915) and X(4350) can give the 
concrete information of their quantum numbers with model 
independent. Additionally, in this work, we also predict the 
decay behaviors of the two remaining second radial excited 
P-wave charmonium states x"n an d x' c \ ■ These findings are 
expected to be revealed in future experiment. 

This study can be extended to include the theoretical study 
of the hidden-charm decay, the radiative decay and the double- 
photon decay of X(3915) and X(4350), which will provide us 
valuable information of their underlying structure. 

Note added: Recently, a theoretical work using QCD sum 
rule lfl9tl shows that it is not possible to describe the X(4350) 
structure as a U + D*Z)* molecular state, which supports our 
effort to explain X(4350) under the conventional charmonium 
to some extent. 
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